Plasmonic excitations in noble metals: The case of Ag. 
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The delicate interplay between plasmonic excitations and interband transitions in noble metals is 
described by means of ab initio calculations and a simple model in which the conduction electron 
plasmon is coupled to the continuum of electron-hole pairs. Band structure effects, specially the 
energy at which the excitation of the d-like bands takes place, determine the existence of a sub- 
threshold plasmonic mode, which manifests itself in Ag as a sharp resonance at 3.8 eV. However, 
such a resonance is not observed in the other noble metals. Here, this different behavior is also 
analyzed and an explanation is provided. 



X 



I. INTRODUCTION 

In a recent letter ,0 Ku and Eguiluz have shown that 
band structure effects are responsible for the positive dis- 
persion of the plasmon line width in K. Such a behavior 
had been -previously found by vom Felde, Sprosser-Prou, 
and Finkn in electron energy loss experiments. They 
interpreted the result as due to short range electron- 
electron correlations beyond the random phase approx- 
imation (RPA). However, Ku and Eguiluz explain the 
result within the RPA in terms of the extra plasmon de- 
cay channels provided by a manifold of unoccupied d- 
symmetry bands. The effect of these localized bands is 
overlooked when treating K as an homogeneous electron 
system in the jellium model. Therefore, band structure 
effects turn out to be decisive in determining, not only 
the existence interband transitions, but also the proper- 
ties of plasmonic excitations. Indeed, ab initio calcula- 
tions of the dynamical response have properly described 
the experimental plasmon dispersion in a variety of sim- 
ple metals.Bcl On the other hand, the jellium model still 
provides a qualitative, though not accurate, description 
of the energy loss spectrum of these systems in the long 
wavelength limit. 

The energy loss spectrum can be obtained either by 
doing electron energy loss spectcpscopya (EELS), inelas- 
tic scattering of X-ray photonsE^^IXSS), or indirectly by 
means of optical measurements. liil In noble metals this 
spectrum has a complicate structure, which bearjidittle 
resemblance to that of simple metals. It is knowr£3 that 
such a structure has to do with the existence of a manifold 
of d-like bands a few eV below the Fermi level. These lead 
to important deviations from even the crudest predic- 
tions for an homogeneous electron liquid. Thus, whereas 
for Ag a narrow resonance is observed at an energy w 3.8 
eV, no well-defined plasmonic excitations appear in the 



spectra of Cu and Au. 

The purpose of this paper is to investigate the structure 
of the energy loss spectrum in the long wavelength limit, 
where plasmonic excitations are expected to be impor- 
tant. In particular, we have focused on Ag for our study, 
but some of our conclusions do also apply to Cu and 
Au, and could be extended to similar metallic compounds 
with fully occupied d-like bands close to the Fermi level. 
The spectrum of_Cu was studied from first principles in 
a previous work.El In the present work, we show that the 
structure of the energy loss spectrum of Ag up to ~ 10 
eV can be understood in terms of a simple model. In 
particular, we are concerned with the role played by the 
relative values of the Drude plasma frequency and the 
threshold energy for the excitation of the d-like bands. 
After clarifying the relevant elements for the existence 
of a plasmonic excitation below the interband excitation 
threshold, numerical experiments in an ab initio frame- 
work for calculating the dielectric matrix of Ag are per- 
formed to confirm this picture. Furthermore, compari- 
son of the ab initio results with experimental data gives 
good agreement for excitations well above threshold, be- 
ing able to reproduce all the structure in the loss spec- 
trum up to w 30 eV. The underestimate of the threshold 
energy in local density calculations also shows up in our 
calculations, giving a lower plasma frequency than the 
experiments (see below for a more detailed discussion). 

The question of how band structure alters the proper- 
ties of plasmons in metals has been previously addressed 
by a number of authors. Wilsono has studied the shift 
in the Drude plasma frequency by the presence of an 
optical absorption band that lies above (below) the free 
electron plasma frequency. He has found that the plasma 
frequency is shifted downwards (upwards) the effect of 
the absorption band. More recently, Sturmlia has studied 
within the RPA the corrections to the dielectric function 
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of an homogeneous electron gas due to a weak lattice 
potential. He has shown that the lattice potential can 
provide the necessary momentum for a plasmon to decay 
into an electron-hole pair below the Landau cutoff.E3 In 
other words, when taking into account the lattice struc- 
ture by folding the free electron bands the first Brillouin 
zone (BZ), a plasmonic excitation can decay by produc- 
ing an interband transition. As a consequence, even in 
th RPA, plasmons acquire a finite lifetime. Core polar- 
ization effects on the plasma frequency have been inves- 
tigated by Zaremba and SturmJia They treated the core 
polarizability as the superposition of the polarizability of 
the isolated lattice ions. The polarization of core elec- 
trons lowers the plasma frequency from its Drude value, 
bringing it closer to the experimental value. 

Recently, the dielectric matrix as well as the energy 
loss spectrum have been evaluated numerically within 
a first principles framework. In these calculations, one 
usually starts from a set of LDAt3 orbitals and en- 
ergies, using time-dependent density functional theory 
(TDDFT)cj an approximation to the density correla- 
tion function can be calculated. Hence, The energy loss 
spectrum is obtained from the imaginary part of the den- 
sity correlation function, as dictated by the fluctuation- 
dissipation theoremEl To mention a few examples of the- 
ses studies, n besides the above mentioned work by Ku 
and Eguiluza, Quong and EguiluzB have thus investigated 
the anisotropy in the plasmon dispersion of aluminumE3. 
Aryasetiawan and Karlssonta have studied the excitations 
in the energy loss spectra of Li. The negative dispersion 
of bulk plasmons in Cs has been explained by Fleszar 
and co-workers.u In general, specific details of the band 
structure of metals as the presence of band gaps in the 
conduction band region are needed to describe the exper- 
imentakenergy loss spectra in some simpleaQ and noble 
metals. u 

The behavior of nobler-metals has been considered by 
Ehrenreich and Phillipp.EI These authors used the data 
obtained in experimental measurements of the reflectiv- 
ity from Cu and AgcJ to obtain the optical and energy 
loss spectra. They separated the Drude and interband 
contributions to the optical response, thus shedding light 
on the reason why the energy of the plasma resonance in 
Ag is shifted down in energy. In a later work, Cooper, 
Ehrenreich and Phillippo extended this analysis to Au. 
They touched upon the question of assigning interband 
transitions to the most relevant features observed in the 
optical specttum. To further clarify this point, Mueller 
and PhillipscJ performed, within the RPA_and based on 
band structure calculations by BurdickJ£3 a numerical 
calculation of the imaginary part of interband contribu- 
tion to the optical response of Cu. Thus, they could give 
a more correct identification of the interband transitions 
that contribute to the optical spectrum at a given en- 
ergy. Their interpretation is consistent with the optical 
data obtained in a number .of reflectivity measurements 
from Ag and Ag alloys E3l3 

As we have mentioned above, the energy loss function 



of Cu has been recently evaluated from first principles 
by Campillo, Ruhio, and Pitarke.B These authors solved 
the Kohn-ShamEj equations on a plane wave basis using 
a norm conserving pseudopotential.tJ Subsequently, they 
evaluated the dielectric and energy loss functions. The 
method is shown to give an accurate description of the 
energy loss spectrum of Cu once the fully occupied d- 
like bands are taken into account as part of the valence 
electron complex. The remaining core electrons do not 
take an active part in the excitations and would only act 
as a global dielectric background as long as the excitation 
energy remains small compared to their binding energy. 
In the present work, we obtain the energy loss spectrum 
for Ag using the same ab initio techniques but, this time, 
we focus on the plasmonic excitations of noble metals. 
In the case of Ag, the appearance of a plasma resonance 
at ss 3.8 eV can be also qualitatively understood with 
help of a simple model, which provides a new perspective 
on how a subthreshold plasmonic mode can appear, and 
what are the parameters controlling the existence of this 
type of excitations. 

The outline of the paper is as follows. In the following 
section, we discuss a model which provides a qualitative 
explanation for some features in the energy loss spectrum 
up to ~ 10 eV. The details of the ab initio computation of 
the spectrum can be found in Sect. [II. A description and 
discussion of the results is given in Sect. IV. Finally, the 
main conclusions of the present paper are summarized in 
Sect. [v|. Atomic units (h — e 2 = m e = 1) are used in all 
mathematical expressions. 



II. A SIMPLE MODEL 

In order to estimate the plasma frequency in Ag the 
classical expression uj p = v '47m can be used. If we set 
n = n c , the density of conduction electrons, i.e., those 
coming from the atomic 5s orbital, we find that u> c v — 8.98 
eV. However, if we furthermore include those coming 
from the atomic 4d orbitals it turns out that u> p = 32.6 
eV. In both cases, we obtain values clear disagreement 
with the experiment, which gives a value of 3.8 eV. Thus, 
as we have remarked in the previous section, Ag does 
not seem to behave as a "simple" metal whose plasmonic 
excitations can approximately understood within the jel- 
lium model. Ehrenreich and Phillippll3 considered this 
issue, concluding that it is also necessary to take into 
account the effect of electrons coming from the atomic 
Ad orbitals in a more subtle way. In fact, it is necessary 
to account for the full band structure of Ag in order to 
explain what is observed in the experiments. In what fol- 
lows, we shall study how the interband transitions affect 
the collective motion of conduction electrons, leading to 
the existence of a subthreshold plasmonic excitation. 

The usual way to find out the energy plasmonic exci- 
tations goes as follows. Disregarding crystal local field 
effects, plasmonic excitations correspond to zeroes in the 
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macroscopic dielectric function. In the |q| = limit we 
have that: 

e(O)=0, (1) 

for same complex fl. Using the Sellmeyer-Drude expres- 
sionllll for c(ui), and neglecting damping effects, 



e ( W ) = i-^§!+* e (*)( W ), 



(2) 



where last term corresponds to the interband contribu- 
tion: 
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U„>0 



The condition in Eq. (hi) can be written as 
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{lu°) 2 +n 2 5e {tb \VL) = 0. 



(3) 



(4) 



The oscillator strength, /„ in Eq. (|J), is a measure of the 
effective number of electrons that participate in the exci- 
tation of a given mode 'n', with energy w„. The dominant 
contribution at frequencies below oj^ comes from the ex- 
citatioB-pf conduction electrons. The Drude plasma fre- 
quency!^ ujp — 9.2 eV, is slightly higher than ujp. The 
ratio of the two plasma frequencies is usually expressed 
in terms of the "optical" mass, m ovt (recall that m e = 1 
in atomic units): 



UJ 



"opt 



0.95 



(5) 



for Ag and m opt — 1.35 for CuJll In the RPA, m opt cor- 
responds to the effective mass of electrons averaged over 
the occupied part of the conduction band. 

Eq. (||) can be read in the following way. We assume 
that the term Q 2 5e^ b ' (f2) plays the role of a "self-energy" 
for the Drude plasmon. In simple metals, this term is 
usually small and varies slowly with to. This leads to a 
small shift of the Drude frequency. The plasma frequency 
also acquires a small imaginary part, which accounts for 
the fact that the plasmon state is degenerate in energy 
with other configurations of the system involving one or 
several electrons excited from one band to another. This 
idea can be easily realized in a model in which the plas- 
mon is viewed as a boson coupled to an energy contin- 
uum. In an independent electron picture, the continuum 
corresponds to the electron-hole pair excitations. A sim- 
ple model Hamiltonian that can be set up to describe this 
situation is the following : 

H = E b l b « 

q 

+ TO L mfc >q 



+ 9k,q (bq + &tq)("4,q + m_ fc) _q). 



(6) 



In this model, plasmons are presented as a field of oscil- 
lators carrying momentum q. Thus, 6 q and 6 q are plas- 
mon creation and annihilation operators, respectively. To 
keep the model as simple as possible, but retaining all the 
important physical effects, the energy continuum is also 
represented by a set of harmonic oscillators whose quanta 
are created and annihilated by the operators m k , TO/^q. 
The label k stands for the additional quantum numbers 
(spin, relative momentum, band indices...) carried by a 
given continuum mode. The operator m_/j._ q (tj^ _ ) 
anhilates (creates) a continuum mode with all the mo- 
mentum quantum numbers reversed. The assumed inter- 
action between plasmons and the energy continuum con- 
serves the momentum q. Terms connecting plasmons and 
modes with momentum differing by a finite reciprocal lat- 
tice vector are therefore neglected. This approximation 
amounts to disregarding crystal local field corrections, 
which in general is a good approximation for metals. In- 
deed, in the frequency range in which we are interested, 



namely uj 



tP 



,s it can be seen from the ab initio 



calculations in CuEl and Ag (see below), the local field 
effects at small q amount to a few percent correction- .— . 

The above model, Eq. (^|), can be solved exactlyB^H 
Indeed, it is related to the models introduced by FanolE 
and AndersonEil to study the coupling of a discrete state 
to a continuum. More precisely, the present model can 
be thought of as an extension of Gadzuk's workEd for a 
localized vibrational mode in an electron liquid. Thus, 
the energies of the plasmonic modes are given by the 
solutions of the following equation: 



n 2 -(LO») 2 + 2^M(q,n) =0, 



(7) 



where M(q, uj) is the plasmon self-energy. It can be ex- 
pressed in terms of its imaginary part: 

Im M(q,w) = |5fc, q | 2 - w*, q ) 

k 

-$(W+W*,q)], (8) 



by means of the expression: 

da/ ImM(q,u/) 



M{q,Lu) 



7T u) — uj — irj 



(9) 



with 77 — s- + . 

Working by analogy with Eq. (ji|), leads us to tenta- 
tively identify 



M(|q| -> 0,w) 



2< 



(10) 



If the RPA expression of Se^ (uj) is used, we can express 
it in terms of single-electron orbitals, |nk), and energies, 



ImM(|q| ->• 0,w) 



2n 2 



E E /(^) 



x [1 - /(e mk )] |(mk|p|nk)| 2 (5(w - e mk + e nk ), (11) 
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where V is the volume of the system and p the momen- 
tum operator. We have assumed cubic symmetry and a 
local pseudopotential; if the pseudopotential is not local, 
the matrix element of p must be replaced by a more com- 
plicated expression. For lo = lo^ this equation matches 
the Golden Rule result for the plasmon line width in the 
|q| = limit,Ea which confirms the identification made 
above in Eq. jl0|). For arbitrary lo, however, it is better 
to look at Im M(q, lo) as a weighted density of the con- 
tinuum modes, cf. Eq. (||). Thus, it follows from Eqs. (0) 
and (^) that if to® lies above (below) the region where 
Im M(q,u)/u)p is large (provided ImM(q,Up)/up be 
small), the plasma frequency will be blue (red) shifted. 
This result is in—agreement with the above mentioned 
work of WilsonO When ImM(q, to®)/ to® is large, we 
will not observe a well-defined plasma resonance in the 
spectrum. This seems to be the case of noble metals due 
to the presence of a manifold of occupied d-like bands, 
which gives rise to a broad band of modes extending from 
the interband threshold energy lot ~ 3.9 eV (for Ag and 
2.1 eV for Cu), to well above the Drude plasma frequency 
(9.2 eV for Ag and 9.3 eV for Cu) . If we assumed that the 
onset of interband transitions in Ag took place so sharply 
as to produce a discontinuity in the weighted density of 
continuum modes: 



damping of the plasmon when the 4d electrons are in- 
cluded in the pseudocore, i.e., they are considered to be 
frozen. In Sect. IV, we shall explain this issue more in 



ImMfjql -> 0,uj) 



(12) 



for lo ~ lot, then ReM(q| — ► 0, lo) would develop a log- 
arithmic behavior ~ — log(w — lot) for lo < lot- This 
means that Eq. ([?]) would have a dampingless solution 
just below the threshold, lot- This polaronic-like so- 
lutionis corresponds to an hybrid plasmonic mode, in 
which the conduction electrons oscillate coherently with- 
out exciting the electrons in the (i-like bands. These are 
just polarized by the electric field set up by plasmonic 
mode, so that the plasma frequency is lowered. This 
point of view|-agrees well with more phenomenological 
approaches ,E3Ea which simply assume that the red shift 
in the plasma frequency can be accounted for by an effec- 
tive dielectric function e^, for the electrons in the (i-like 
bands so that lo* = to® / ' ^ftd ~ 3.8 eV. 

Additionally, there is another complex solution of 
Eq. (0), for which the conduction electrons move as an 
overdamped oscillator, rapidly decaying into an electron 
hole-pair. For |q| — + 0, this excitation corresponds to an 
optical interband transition. Two types..of optical tran- 
sitions may occur above the threshold. Ej Two examples 
have been indicated by arrows in Fig. |l|, which represents 
the calculated first-principles LDA band structure of Ag. 
Type A corresponds to the excitation of an electron from 
one flat d-like band to an unoccupied state in the conduc- 
tion band or above. The other type (B), whose threshold 
occurs at approximately the same energy in the three 
noble metals, corresponds to the excitation from an oc- 
cupied state in the conduction band (near L point in the 
BZ) to the flat part of an unoccupied band or above. 
This type of transitions seems to be responsible for the 



detail. These two types of transitions are important be- 
cause they often involve flat bands, which may lead to 
van Hove singularities in the joint density of states. 

In the case of Ag, the onsets for the transitions of 
type A and B nearly overlap in energy. This pro- 
duces a sharp interband onset and gives rise to the pro- 
nounced peak that line (a;) exhibits around 4 eV. From 
Eq. @, it follows then that ImM(|q| -> 0,to) will 
also have a maximum at the same energy. The onset 
of the interband transitions at lot is not as sharp as sug- 
gested by Eq. (|l2|). However, it still leads to a zero of 
n 2 - (lo°) 2 + 2t^ReM(|q| -► 0,O) for a real ft, which 
occurs below the threshold, where the density of contin- 
uum modes is very small. This point is illustrated by 
Fig. ||, which displays the graphical solution of Eq. (0). 
This provides an explanation for the existence of a sub- 
threshold plasma resonance in Ag. However, it remains 
to be explained why a similar phenomenon is not ob- 
served in the spectra of Cu or Au. The difference stems 
from the specific details of the band structure and, in 
particular, from the value of the threshold energy. In 
Fig. |[ we have plotted the graphical solution of Eq. (0)in 
the case of Cu. The threshold for the excitation of (i-like 
bands occurs this time at a lower energy (2.1 eV) because 
the onset for the interband transitions of type A and B 
does not overlap. Comparison of Figs. ^ and ^ shows 
the important role played by the value of the thresh- 
old energy, cf. Eq. (|l0|), in enhancing the maximum of 
ImM(|q| -> 0,lo), and thus in ReM(|q| -> 0,lo). There- 
fore, for Cu no zero of Eq. (^) close to the real axis exists 
below the threshold, where this plasma mode would have 
a small density of continuum modes to decay. A similar 
conclusion can be drawn for the case of Au. 

To sum up, band structure effects modify the original 
Drude plasmon, which turns out to be no longer a well 
defined plasmonic excitation. Thus, it shows up in the 
Ag spectrum as broad peak with a maximum around 8 
eV. Moreover, a combination of a high value for lot along 
with the sudden onset of interband transitions produces a 
narrow plasma resonance in Ag. In this respect, our cop, 
elusions agree with those of Ehrenreich and Phillipp.E3 
Furthermore, the line width of this "delicate" feature in 
the energy loss spectrum also depends on the sharpness 
of the interband onset. The sharper onset the narrower 
the resonance, which in any case has very little oscillator 
strength. This makes this plasmonic excitation very sen- 
sitive to changes in the band structure, impurities or de- 
fects. This is consist(p±^ith the experimental data from 
alloying experimentsJlJ'oEj which indicate that the plas- 
mon is strongly damped by increasing the concentration 
of the other component in the alloy. Finally, it is worth 
pointing out that the existence of this resonance also de- 
pends on the value of the Drude plasma frequency. If the 
value of LOp is increased, for example by increasing the 
electron density in the system, the resonance could dis- 
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appear from the spectrum. This prediction is confirmed 
by our ab initio results (shown below) for the energy 
loss spectrum of Ag under externally applied hydrostatic 
pressure. 



III. DETAILS OF THE AB INITIO 
CALCULATIONS 

By measuring the energy lost by electrons or X-ray 
photons in their interaction with matter, one can ac- 
cess the dynamical structure factor ,-S(q, lo) of a system. 
The fluctuation-dissipation theoremcil allows us to relate 
5(q, uj) to the density correlation function: 

X(r,r',t) = -i6(t) ([5n H (r,t),Sn H (r',0)]), (13) 

where 5riH(r,t) = n#(r, t) — (n#(r, £)), and suffix H 
means that operators are in the Heisenberg pictured; the 
brackets mean that the average is taken over the ground 
state of the system. Now, time-dependent density func- 
tional theorycd (TDDFT) provides us with a method to 
compute this correlation function by solving the following 
integral equation: 



XGG'(q,w) = XGG'(^ ( 



E 

Gi .G2 



[w Gl G 2 (q) +^c 1 G 2 (q, w )] XG 2 G'(q,w). 



(14) 



Here, we have exploited crystal symmetry and time trans- 
lation invariance by introducing: 



XGG'(q,w) = ^ J dt d 3 r d 



3 r' e iut e ~ i (i+ G )' r 



xe 



i(q+G')-r' 



X(r,r',t), (15) 



where q belongs to the BZ, and G, G' to the reciprocal 
lattice. Similar definitions hold for the other functions 
that appear in Eq. (|l4|). For example, 



fGG'(q) = 



An 



G 



-JOG.G' 



(16) 



corresponds to the Fourier components of the Coulomb 
interaction between metal electrons. The function 



XGc(q^) 



kGBZ n,m 



mk+q ) 



■ £nk - £mk+q + 



x(nk|e- l(q+G) - r |mk + q) 
x(mk + q|e i ( q+G ')- r |?^k) (17) 

with 77 — ^ + , is the density correlation function for a 
system of independent particles, with single-particle en- 
ergies e„k and orbitals |nk). These are eigenvalues and 
eigenstates of a single-particle Hamiltonian. The poten- 
tial epergy in this Hamiltonian is given by the Kohn- 
Shamc3 potential. In our case, the JCohn-Sham potential 
has been calculated within LDALLSEij. 



We have used a plane wave basis to expand the Bloch 
eigenstates. For this procedure to be computationally 
efficient, it is necessary to replace the lattice potential 
by a superposition of ionic pseudopotentials having the 
same scattering properties. However, as electrons in the 
outer atomic Ad shell play an important role in the re- 
sponse of noble metals, they are treated as a part of the 
valence electrons. The remaining core electrons are re- 
placed by a frozen pseudacore using a norm-conserving 
pseudopotential scheme£j Therefore, electrons included 
in the pseudocore will not respond to an external per- 
turbation. This is not an important drawback as long as 
the excitation energy remains below the typical excita- 
tion energies of core electrons ~ 100 eV. 

Including Ad electrons in the valence produces a pseu- 
dopotential that requires a relatively high number of 
plane waves to get good convergence. Here, we have used 
a scalar relativistic pseudopotential generated according 
to the scheme of Troullier and MartinsEj It requires 
« 1, 400 plane waves per Bloch state up to a cutoff energy 
of sa 1.1 keV in order to obtain a well converged ground 
state energy (less than 0.01 eV), and single-electron en- 
ergies and orbitals up to « 100 eV. The lattice constant 
has been set to the experimental value of the FCC crystal 
structure of Ag, namely a = 4.09 A. We have also solved 
the Kohn-Sham equations with a = 3.70 A, which corre- 
sponds to an externally applied hydrostatic pressure of 
63 GPa. 

The RPA is equivalent to solving Eq. (0) with K xc = 
0. Therefore, all the exchange correlation effects beyond 
the RPA are contained in the kernel K G c G ,(q, ui). Con- 
sistently with the use of LDA eigenstates and eigenvalues 
in the calculation of x > this function is approximated by 
the following expression: 



K xc 

A GG' 



1 

V 



d r e 



-i(G-G')-r 



/xc(r), 



where 



/*c(r) 



dV xc (n) 



dn 



(18) 



(19) 



»(r) 



Here V xc is the exchange correlation part of the Kohn- 
Sham LDA potential, and n gs (r) is the ground state den- 
sity. This approximation is usually called adiabatic local 
density approximation (ALDA). However, it is observed 
that for metals in the long wavelength limit, ALDA and 
RPA yield veW|-similar results. n This has been found in 
simple metalsjd'cl and also in GuB Thus, when computing 
the spectrum, we have used the RPA. 

Instead of solving Eq. ( p"7j ) , we have evaluated the RPA 
dielectric matrix: 

£ gg' (q, = 5 G ^ G , - v GG > (q)xcc' (q. w )- ( 20 ) 

In the RPA, the dynamic structure factor 5(q, u>) is pro- 
portional to the loss function, Im [— e^, (}(q, w)] , which is 
the G = G' = element of the inverse dielectric matrix. 
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To calculate the dielectric matrix we have used a cutoff of 
41 eV in G, G', which amounts to using 59 x 59 matrices. 

To evaluate x° using Eq. ( |T7j ) , we have cut off the sum 
over n and m using bands up to an energy of ~ 100 eV. 
The integration over the first BZ has been_performed 
using a 20 x 20 x 20 Monkhorst-Pack mesho in the ir- 
reducible BZ, which allows to deal effectively with the 
topology of the Fermi surface of Ag. Indeed, this is a very 
dense mesh and the calculation of \° becomes certainly 
time and memory consuming. However, if one is to re- 
solve the very narrow plasma resonance that Ag presents 
near the interband excitation threshold, it is then neces- 
sary to work with such a fine mesh. 

In order to further reduce the numerical damping r\ 
we have also computed \° on t ne imaginary frequency 
axis by making the replacement lo + irj — » iv in Eq. (17). 
Analytic continuation to the real axis has been carried 
out by using a Pade approximant.c£l This procedure al- 
lows us to use a numerical damping as small as 0.001 
eV. This is to be compared with the 0.2 eV value used 
in the real frequency calculations also presented in this 
work. However, we shall show the results of boht types 
of calculations together, to indicate that one must be 
very careful with the analytic continuation procedure us- 
ing Pade approximants. Thus, the results obtained using 
this method tend to be smoother than what is obtained 
with a calculation with real frequencies and larger values 
of rj. Thus, it usually happens that some features of the 
spectrum, as obtained in a real frequency calculation, are 
completely missing from the analytically continued result 
(see following section). In fact, although we have used 
Lentz algorithmcil to evaluate the corresponding contin- 
ued fraction, we have been unable to ensure that the 
procedure converges by increasing the number of points 
over the same interval of v. When fitting the Pade ap- 
proximantE3 to the calculated values of ec,G(q, iv), we 
have found that the results are unstable with respect to 
increasing the number of points. Thus, we have obtained 
differences of the order of 10 % when the number of points 
over the same interval was simply doubled. Finally, we 
have decided to use the modulus of the remainder of the 
truncated continued fraction to estimate the error, and 
always compare the outcome with a real frequency cal- 
culation prior to consider it as reliable. 



IV. RESULTS 

In Fig. [I] we have plotted the band structure of Ag as 
calculated within the LDA. Two examples of the main 
types of interband transitions (A and B) which couple to 
plasmons at energies u> ~ lo® are indicated by arrows. 
In the LDA, we obtain a threshold energy for type A 
transitions w 3 eV, whereas experimentally it is found 4 
eV. Since LDA eigenvalues enter the expression for x j 
this turns out to be the reason for underestimating the 
threshold frequency, lot, as shown in Figs. || and [|. An 



underestimate of this kind was already observed for Cu, 
where lot ~ 1-5 eV in contrast with the experimental 
value of 2.1 eV. The source of such an underestimate 
may be cither a failure of the LDA or an indication that 
a frequency dependent and perhaps nonlocal K xc may 
be required to fully account for the experimental value 
of lot- In our opinion, however, it is related to a failure 
of the LDA, which may not be taking into account the 
strong correlations occurring amongst the electrons in 
the flat d-like bands of Ag. Therefore, a more detailed 
description p£ correlation for the d-like bands seems to 
be required a 

Fig. [I] presents our results for the energy loss func- 
tion for q = (0.05 0)(27r/a) as obtained from the LDA 
band structure. Since our interest is focused on the plas- 
monic excitations, which occur for small q, we need not 
consider exactly the |q| = limit. Thus, for the small 
q considered, the theory developed in Sect. || still ap- 
plies. The continuous line in this figure corresponds to 
Im — £oo(q, lo)] , whereas the dashed line corresponds to 
Im — l/eoo(q 5 w)], i.e. with and without including crys- 
tal local field effects respectively. As we have already said 
in Sect. 0, crystal local field effects are at most a cor- 
rection ~ 10 % in the considered frequency range. Dots 
are the experimental dataEl obtained by means of optical 
measurements (i.e., for |q| — ► 0), and should be compared 
with the dashed line. As a consequence of the LDA un- 
derestimate of lot, the plasmon peak appears at-a lower 
energy w 3.4 eV than observed experimentallyE3 ss 3.8 
eV. This peak also looks more damped either because 
the high value of the numerical damping employed in the 
calculation and because of the LDA underestimate of lot- 

From Figs. [| and || we can conclude that: (i) The 
Drude plasmon becomes so broadened by its coupling 
to interband transitions that we can no longer speak of 
a well denned excitation. The high energy tail of this 
broad feature around 8 eV is more strongly affected by 
crystal local field effects. This seems to indicate that 
those decay processes involve the excitation of an elec- 
tron in a more localized state, probably lying deeper in 
the d-like band. The important role played by the d- 
like bands in destroying the Drude plasmon can be seen 
in Fig. [| where the dashed line corresponds to a cal- 
culation of the energy loss function putting the ten 4<i 
electrons of atomic Ag in the frozen pseudocore. In this 
case, the Landau damping of the Drude plasmon is only 
due to transitions of type B, which have a much smaller 
spectral weight than those of type A. Furthermore, the 
plasma frequency is pushed upwards by this continuum 
of modes, which mainly lies below lo® . However, when 
electrons from the d-like bands are considered as active, 
the effect goes in the opposite direction thus shifting the 
maximum of the Drude peak from lo® eV to ~ 7.5 eV; 
(ii) In agreement with the experiments, a double hump 
structure in the region between 15 to 30 eV is found. This 
is a band structure effect, which tur-ns out to be a general 
feature in the noble metal systems.! As it can be seen in 
Fig. [| although the calculation using a Pade approxi- 
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mant (dot-dashed line) gives a reliable description of the 
low frequency part, it fails to reproduce the double hump 
structure. This seems to be a rather general feature of 
analytic continuation by Pade approximants. Whenever 
the function to be continued presents two poles close to 
each other and having a large imaginary part, the values 
for the function on the real axis correspond to an aver- 
age of the two poles, and a single broad peak shows up 
in the outcome of the analytical continuation, (iii) Al- 
together, the ab initio LDA calculations presented here 
give a good description of the high energy features (i.e. 
above 10®) of the energy loss spectrum of Ag. These in- 
clude the Landau damping of the Drudc plasmon and the 
double hump structure mentioned above. This could be 
expected on the basis that the spectral weight in the high 
energy region of the spectrum comes, to a large extent, 
from the excitation of electrons to bands well above the 
Fermi level, which are expected to be more accurately re- 
produced by the LDA, as it is an approximation suitable 
for an homogeneous electron system. However, the LDA 
underestimate of the threshold energy is responsible for 
failing to reproduce the low energy structure. This issue 
is further investigated below. 

In Fig. |6|, we show the results for the real and imag- 
inary parts of eoo(q, lo) for q = (0.05 0)(27r/a). The 
same remarks made above apply in this case as well. 
The experimental results, full circles for the imaginary 
part and open squares for the real part are from Ref. |43|. 
In this figure we can see more easily that the onset for 
the excitation of interband transitions occurs, not only 
at a lower energy, but that it is less sharp as well. As it 
has been remarked in the previous section, this leads to 
a more damped plasmonic mode below the threshold. 

In order to investigate the influence of the interband 
threshold energy on the existence of a subthreshold 
plasma resonance, we have performed a series of numeri- 
cal "experiments" . The motivation of these experiments 
is to understand why Ag presents a narrow plasma res- 
onance in its energy loss spectrum, while Cu (and Au) 
does not, despite the similarity of their band structures. 
To modify the threshold energy lot without opening un- 
physical gaps in the band structure, we have scaled by 
a factor a the separation between the Fermi energy (sf) 
and the LDA eigenvalues of the occupied states: 

£ nk = £ f + a(e nk - e F ). (21) 

Thus, using this new set of eigenvalues, we have calcu- 
lated x° ancl hence obtained the energy loss function. 
This rather ad hoc procedure does not modify the wave 
functions accordingly. Furthermore, it introduces a dis- 
continuity of Vk£^ k at the Fermi surface, so that it might 
lead unphysical results. However, as shown in Fig. [?], this 
seems not to be the case, at least for lo ~ lot and q fi- 
nite. As predicted in Sect, [n], the subthreshold plasmon 
peak is enhanced and shifts up in energy as lot increases. 
In this figure, continuous lines correspond to calculations 
for real frequencies with n = 0.2 eV, while dashed lines 



correspond to the analytic continuation using a Pade ap- 
proximant from the results obtained on the imaginary fre- 
quency axis. Both calculations show the same tendency. 
On the other hand, when the threshold energy is fur- 
ther decreased, from its LDA value ~ 3 eV, to lot ~ 1.5 
eV, the plasma resonance disappears and the spectrum 
becomes more "Cu-like" . 

Finally, in Fig. || we show the energy loss spectrum of 
Ag under an external pressure of 63 GPa. As in previous 
cases, the continuous line corresponds to a calculation 
using real frequencies while the dashed line is the result 
of an analytical continuation using a Pade approximant. 
Notice that in this case the lattice constant is a press — 
3.7 A so that now q = (0.05 0)(27r/a press ), corresponds 
to a slightly larger momentum (|q pre ss| = 8.5 x 10~ 2 A 
vs |q| = 7.7 x 10~ 2 A" 1 ). 

Decreasing the lattice constant changes the width of 
the bands so that they disperse more rapidly with k (i.e., 
the band width becomes larger). Furthermore, electronic 
density also increases. This shifts the maximum previ- 
ously at Lo max m 7.5 eV to w™%* a » 8.8 eV. Since this 
maximum is related to the Drude plasmon, whose bare 
frequency scales with the lattice parameter as a~ 3 / 2 , we 
could understand the shift (neglecting the change in the 
optical mass) as due to the increase in the electronic den- 
sity. Indeed, this argument works reasonablely well as 
(a press /a)' 3 / 2 = 1.16, and w™ a e x ss / 'Lo max m 1.17 . More- 
over, as the bands become wider, specially those corre- 
sponding to excited states, more oscillator strength is 
transferred to higher energies. This decreases the den- 
sity of modes at the Drude frequency, and thus the cor- 
responding peak is now less broad. What is more, in this 
spectrum the subthreshold plasmon is not present any 
more. Again, this is due to the increased electronic den- 
sity, which yields a higher value for the Drude plasma 
frequency so that Eq. (0) has no longer a solution close 
to the real axis. The behavior under pressure predicted 
here remains to be confirmed by experiments. 

V. CONCLUSIONS 

Here, we summarize the main conclusions of the 
present work. We have shown that a first principles cal- 
culation of the energy loss spectrum is able to reproduce, 
within the RPA, the high energy structure in the spec- 
trum of Ag. However, since the LDA underestimates the 
interband energy threshold, the agreement is not as good 
as for the high energy region around 4 eV. For example, 
the plasma resonance observed experimentally at 3.8 eV 
appears as a more damped peak at 3.4 eV. In our opin- 
ion, this underestimate may be related to a failure of 
the LDA to account for electron correlation effects that 
take place in the flat d-like bands of Ag. In fact, an ad 
hoc scaling of the occupied bands so that the threshold 
energy is brought closer to the experimental value gives 
rise to a better defined subthreshold plasmonic excita- 
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tion. We have also shown the important role played by 
the d symmetry bands in determining the existence of the 
plasma resonance at 3.8 eV, just below the interband ex- 
citation threshold. By using a simple model, it is shown 
that the appearance of this excitation has to do with the 
sharpness of the interband onset and the relative values 
of the threshold energy and the Drude plasma frequency 
Since it contains very little oscillator strength, it is very 
sensitive to impurities, defects, and changes in the band 
structure. Thus, we have demonstrated by performing 
numerical experiments that its presence depends on the 
value of the threshold energy. Finally, we have also com- 
puted the energy loss spectrum of Ag under pressure, and 
found that its features can be aslo understood in terms 
of the simple model presented here. We remark that 
the techniques used and analysis carried out in this work 
are not specific of Ag, and could be readily extended to 
other metallic compounds with fully occupied d symme- 
try bands close to the Fermi level. 
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FIG. 1. Calculated LDA band structure of Ag. Arrows 
indicate the two examples (A and B) of the types of inter- 
band transitions that couple to the plasmonic excitations in 
the long wavelength limit, for energies above the interband ex- 
citation threshold. The crystal structure is FCC with lattice 
parameter equal to a = 4.09 A. The zero energy corresponds 
to the Fermi level. 

FIG. 2. Graphical solution of Eq. (Q) for the subthresh- 
old plasmonic modes in Ag. Data for 5e^ '(uj) (i.e., for 
Re M(|q| — > 0, w)) are from Ref. [H| The continuous line 
cuts the dashed line at the point indicated by the black dot. 
This point lies just below the threshold, where the density 
of continuum modes is small. As a consequence, a narrow 
plasma resonance appears in the spectrum. 



FIG. 5. Energy loss spectrum of Ag 

for q = (0.05 0)(27r/a), |q| = 7.7 x 10 _2 A"\ The con- 
tinuous line is a calculation for real frequencies with 77 = 0.2 
eV, cf. Eq. dTij). The dot-dashed line, however, is obtained 
by analytically continuing to the real axis an imaginary fre- 
quency calculation, using a Pade approximant. As it can be 
seen, the results obtained by analytical continuation cannot 
reproduced all the features in the spectrum obtained using a 
real lo (but higher numerical damping rj). The dashed line is 
a calculation of the loss function including the Ad electrons 
of Ag in the frozen pseudocore that replaces the true ionic 
core. Then, these electrons cannot be excited, thus showing 
the important role that they actually play in modifying the 
spectrum. Notice that the intensity of this peak is five times 
what is represented here. 

FIG. 6. Dielectric function eoo(q, w) of Ag for 

q = (0.05 0)(27r/a),|q| = 7.7 x 10~ 2 A -1 . The continuous 
fine corresponds to the imaginary part as obtained by analyt- 
ically continuing an imaginary frequency calculation, using a 
Pade approximant. The dashed line is the real part obtained 
by the same method. We checked that the results are con- 
sistent with a real frequency calculation (not showed here for 
clarity). Full circles (imaginary part) and open squares (real 
part) correspond to experimental data form Ref. f43J. 



FIG. 3. Graphical solution of Eq. (Q) for the subthresh- 
old plasmonic modes in Cu. Data for <5e' t6 '(u) (i.e., for 
Re M(|q| — » 0, w)) are from Ref. |l^. The continuous line 
does not cut the dashed line, indicating that no subthreshold 
plasmonic modes exist in Cu. A similar situation is expected 
to hold for Au. 

FIG. 4. Energy loss spectrum 

of Ag for q = (0.05 0)(27r/a), |q| = 7.7 x 10" 2 A -1 . The 
continuous line is imaginary part of G = G' = element 
of the inverse dielectric matrix, namely the energy loss func- 
tion Im [— £qq (q, a;)] , whereas the dashed line corresponds to 
(the imaginary part of) one over the G = G' = element of 
the dielectric matrix Im [— l/eoo(q, &>)]. Differences between 
them are due to crystal local field effects. Dots are the ex- 
perimental data from Ref. ||. As they were obtained from 
optical data and do not include crystal local field corrections, 
they must be compared to the dashed line. 



FIG. 7. Results of the scaling procedure described in 
Sect. |[y| The continuous line are the results as calculated 
for real cu with r\ = 0.2 eV, while the dashed lines are ob- 
tained by analytical continuation using Pade approximants. 
When the occupied bands are scaled so that the interband 
threshold energy increases, the subthreshold plasmon appears 
more defined. However, when threshold energy is decreased 
the plasmon disappears from the energy loss spectrum. 

FIG. 8. Energy loss spectrum 

of Ag for q = (0.05 0)(27r/a press ) (|q| = 8.5 x 10" 2 A -1 ) 
with dpress ~ 3.7 A, which corresponds to an externally ap- 
plied hydrostatic pressure of 63 GPa. The continuous line is a 
calculation for real frequencies with 77 = 0.2 eV, cf. Eq. (p"7j). 
The dashed line, however, is obtained by analytically continu- 
ing to the real axis an imaginary frequency calculation, using 
a Pade approximant. 
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